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Abstract  Decisions  about  recreational  beach  clo¬ 
sures  would  be  enhanced  if  better  estimates  of 
surfzone  contaminant  transport  and  dilution  were 
available.  In  situ  methods  for  measuring  fluo¬ 
rescent  Rhodamine  WT  dye  tracer  in  the  surf¬ 
zone  are  presented,  increasing  the  temporal  and 
spatial  resolution  over  previous  surfzone  tech¬ 
niques.  Bubbles  and  sand  suspended  by  breaking 
waves  in  the  surfzone  interfere  with  in  situ  opti¬ 
cal  fluorometer  dye  measurements,  increasing  the 
lower  bound  for  dye  detection  1  ppb)  and  re¬ 
ducing  (quenching)  measured  dye  concentrations. 
Simultaneous  turbidity  measurements  are  used 
to  estimate  the  level  of  bubble  and  sand  inter¬ 
ference  and  correct  dye  estimates.  After  correc¬ 
tion,  root-mean-square  dye  concentration  errors 
are  estimated  to  be  <  5%  of  dye  concentration 
magnitude,  thus  demonstrating  the  viability  of  in 
situ  surfzone  fluorescent  dye  measurements.  The 
surfzone  techniques  developed  here  may  be  ap¬ 
plicable  to  other  environments  with  high  bubble 
and  sand  concentrations  (e.g.,  cascading  rivers 
and  streams). 


D.  B.  Clark  (ISI)  •  F.  Feddersen  • 

M.  M.  Omand  •  R.  T.  Guza 

Scripps  Institution  of  Oceanography,  UCSD, 

9500  Gilman  Dr.,  mail  code  0209,  La  Jolla,  CA 

92093-0209,  USA 

e-mail:  clark@coast.ucsd.edu 

Published  online:  16  April  2009 


Keywords  Surfzone  •  Mixing  •  Dispersion  • 
Fluorometer  •  Turbidity  •  Fluorescent  dye 
measurement 


1  Introduction 

Surfzone  dispersion  is  important  to  many  biolog¬ 
ical  and  physical  processes,  including  the  dilu¬ 
tion  of  contaminated  runoff  (e.g.,  Boehm  2003; 
Grant  et  al.  2005),  phytoplankton  transport  (e.g., 
Campbell  and  Bate  1988),  and  larval  recruit¬ 
ment  (e.g.,  Denny  and  Shibata  1989).  However, 
few  surfzone  tracer  dispersion  studies  have  been 
conducted,  and  measurement  techniques  are  less 
advanced  than  those  used  further  offshore.  As  a 
result,  surfzone  tracer  dispersion  is  poorly  under¬ 
stood.  Increased  knowledge  of  transport  and  dilu¬ 
tion  in  the  surfzone  would  enable  beach  managers 
to  better  minimize  public  exposure  to  shoreline 
pollution. 

Fluorescent  dye  (e.g.,  Rhodamine  WT,  or  flu¬ 
orescein),  optically  measured  with  a  fluorometer, 
is  commonly  used  to  study  tracer  (i.e.,  pollutant) 
dispersion  in  both  marine  (e.g.,  Okubo  1971)  and 
freshwater  (e.g.,  Csanady  1963)  environments. 
Dye  measurements  are  used  to  estimate  eddy 
diffusivities  and  scale  dependence  (e.g.,  Okubo 
1971;  Murthy  1976;  Fong  and  Stacey  2003),  tracer 
transport  (e.g.,  Houghton  1997;  Tilburg  et  al. 
2007),  and  other  quantities.  The  observations  of 
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diffusion  and  mixing  can  be  used  to  calibrate  nu¬ 
merical  models. 

Ideally,  accepted  in  situ  techniques  using  flow¬ 
through  (e.g.,  Pritchard  and  Carpenter  1960)  or 
open-face  (e.g.,  Ledwell  et  al.  2004)  fluorome- 
ters  could  be  adapted  for  surfzone  use.  However, 
surfzone  waves  and  strong  currents  can  damage 
instrumentation,  hamper  deployment,  and  render 
many  oceanographic  instruments  and  platforms 
(i.e.,  boats)  unusable.  Breaking  waves  intermit¬ 
tently  suspend  optically  interfering  bubbles  (e.g., 
Deane  1997)  and  sediment  (e.g.,  Brenninkmeyer 
1976;  Yu  et  al.  1993),  creating  highly  variable 
turbidity. 

Measured  dye  fluorescence  is  altered  by 
many  factors,  including  background  fluorescence 
(Pritchard  1979),  temperature,  and  suspended 
sediment  (Smart  and  Laidlaw  1977).  Suspended 
sediment  increases  apparent  background  fluo¬ 
rescence  without  dye  in  the  water  and  reduces 
(quenches)  concentration  estimates  when  dye  is 
present.  However,  the  effect  of  suspended  sed¬ 
iment  on  dye  estimates  varies  with  sediment 
type,  color,  and  concentration  (Smart  and  Laidlaw 
1977),  and  the  effect  of  suspended  sand  in  the 
surfzone  is  unknown.  The  effect  of  surfzone  bub¬ 
bles  on  in  situ  dye  measurements  is  also  unknown. 
Dye  mass  can  also  be  lost  to  photochemical  decay 
(e.g.,  Suijlen  and  Buyse  1994),  adsorption  (e.g., 
Talbot  and  Boon  1975),  and  other  factors  (e.g.. 
Smart  and  Laidlaw  1977),  but  these  processes  are 
not  expected  to  be  significant  over  the  few-hour 
duration  of  most  surfzone  dispersion  studies. 

Previous  surfzone  dye  measurements  have 
been  limited  to  bottle  samples  taken  either  at 
the  shoreline  or  near  the  visually  estimated  cen¬ 
ter  of  a  spreading  dye  patch  (Harris  et  al.  1963; 
Inman  et  al.  1971;  Clarke  et  al.  2007).  While  sim¬ 
ple  and  inexpensive,  the  number  and  distribution 
of  samples  has  been  severely  limited,  with  maxi¬ 
mum  rates  up  to  two  bottle  samples  per  minute 
(Clarke  et  al.  2007),  and  a  maximum  of  four  si¬ 
multaneous  locations  over  several  hours  (Harris 
et  al.  1963).  These  studies  provided  important 
initial  estimates  of  surfzone  dispersion  but  could 
not  resolve  cross-shore  dye  structure  and  used  few 
(or  single)  dye  patch  or  plume  realizations.  In 
situ  fluorometers  allow  for  higher-frequency  dye 
measurements,  but  have  not  been  implemented 


in  the  surfzone,  possibly  because  of  concerns  that 
suspended  bubbles  and  sand  would  degrade  in¬ 
strument  performance. 

Instruments  suitable  for  in  situ  surfzone  mea¬ 
surement  of  fluorescent  Rhodamine  WT  dye  are 
discussed  in  Section  2,  along  with  possible  mech¬ 
anisms  causing  high  errors  in  some  instruments. 
The  specific  instrumentation  used  here  (Section  3) 
and  testing  procedures  (Section  4)  are  described. 
The  effects  of  surfzone  bubbles  and  sand  on  in  situ 
Rhodamine  WT  measurements  are  characterized, 
and  results  are  given  in  Section  5.  Errors  in  in 
situ  surfzone  fluorescence  are  estimated,  recom¬ 
mendations  are  made  for  error  reduction,  and 
example  field  applications  are  given  in  Section  6. 
Details  of  mixing  and  delay  time  within  the  mobile 
flow-through  fluorometer  system  and  estimates 
of  effective  spatial  resolution  are  given  in  the 
Appendix. 

2  Considerations  for  in  situ  Surfzone  Rhodamine 
WT  Measurement 

Many  commercial  fluorometers  are  available.  Two 
widely  used,  open-face-type  fluorometers  (WET 
Labs  ECO  Triplet  and  Turner  SCUFA)  that  mea¬ 
sure  Rhodamine  WT  dye  fluorescence  D  and 
turbidity  r  were  tested  side-by-side  in  a  natural 
surfzone  without  dye  in  the  water  to  determine 
the  instruments’  responses  to  suspended  bub¬ 
bles  and  sand  from  breaking  waves.  Light  emit¬ 
ted  near  the  dye  excitation  peak  is  fluoresced 
at  a  slightly  longer  wavelength  (e.g.,  Guibault 
1990).  The  intensity  of  the  fluoresced  light  is 
detected  at  an  angle  to  the  excitation  beam  (scat¬ 
tering  angle)  and  used  to  determine  dye  con¬ 
centration  (Fig.  1).  The  WET  Labs  ECO  Triplet 
measures  D  (0.1-400  ppb  range)  with  530  nm 
[30  nm  full  width  half  max  (FWHM)]  wavelength 
excitation,  and  570  nm  fluorescence  detection 
(35  nm  FWHM),  and  measures  r  (0.03-100  ntu 
range)  at  660  nm  wavelength.  ECO  Triplet  D 
and  T  both  use  117°  scattering  angles.  The  Turner 
SCUFA  measures  D  (0.04-400  ppb  range)  with 
530  nm  (40  nm  FWHM)  wavelength  excitation, 
and  600  nm  (40  nm  FWHM)  fluorescence  de¬ 
tection,  and  measures  r  (0.1-400  ntu*  range)  at 
530  nm.  SCUFA  D  and  r  both  use  90°  scattering 
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Fig.  1  Schematic  of  fluorometer  optics  (ECO  Triplet  ex¬ 
ample).  Light  emitted  near  the  excitation  peak  of  Rho- 
damine  WT  is  fluoresced  by  the  dye  at  a  longer  wavelength 
and  measured  by  the  detector.  When  dye  concentrations 
are  low,  bubbles  and  sand  in  the  sample  volume  {center 
shaded  region)  scatter  excitation  light  from  the  instrument 
towards  the  dye  fluorescence  detector.  Overlap  between 
the  instrument  excitation  spectrum  and  the  fluorescence 
detection  spectrum  allows  backscattered  excitation  light  to 
be  falsely  interpreted  as  an  elevated  dye  concentration  in 
the  sample  volume.  When  dye  concentrations  are  not  low, 
bubbles  and  sand  reduce  the  amount  of  fluoresced  light 
reaching  the  detector,  and  measured  dye  concentrations 
are  reduced.  Bubbles  and  sand  also  backscatter  light  from 
the  turbidity  light  source  into  the  turbidity  detector  (not 
shown),  increasing  measured  r 


angles.  ECO  Triplet  and  SCUFA  r,  with  units  ntu 
and  ntu*,  respectively,  are  similar  but  not  equiva¬ 
lent  due  to  instrumental  differences  (Section  3.2). 

The  instruments  were  mounted  1  m  apart,  so 
they  were  close  enough  to  measure  similar  tur¬ 
bidities  but  far  enough  apart  to  eliminate  inter¬ 
ference.  Both  the  ECO  Triplet  (Fig.  2a)  and  the 
SCUFA  (Fig.  2b)  recorded  spurious  dye  mea¬ 
surements  (nonzero  dye  measurements  when  no 
dye  is  present)  in  response  to  r  from  surfzone 
bubbles  and  sand.  The  ECO  Triplet  recorded  very 
little  spurious  dye  (<  0.8  ppb),  while  the  SCUFA 
recorded  large  spurious  dye  measurements  (up 
to  17  ppb)  that  were  highly  correlated  with 
turbidity  r. 

The  SCUFA  and  ECO  Triplet  measure  Rho- 
damine  WT  fluorescence  with  different  opti¬ 
cal  elements  and  angles.  Ideally,  narrow-banded. 


nonoverlapping  optical  filters  on  the  source  and 
detection  elements  prevent  excitation  light  from 
being  detected  as  dye.  However,  nonideal  opti¬ 
cal  filters  allow  excitation  light,  backscattered  off 
bubbles  and  sand  (Fig.  1),  to  pass  through  the 
detection  filter  and  to  be  measured  as  elevated 
dye  concentrations.  The  SCUFA  optical  filters 
have  a  slightly  wider  pass-band  than  the  ECO 
Triplet,  possibly  increasing  SCUFA  spurious  dye 
measurements.  However,  the  SCUFA  filter  cen¬ 
ters  are  also  10  nm  further  apart  than  the  ECO 
Triplet,  slightly  reducing  this  effect.  The  intensity 
of  light  scattered  by  surfzone  bubbles  and  sand  in¬ 
creases  with  a  reduction  in  scattering  angle  (Zege 
et  al.  2006).  Therefore,  the  smaller  SCUFA  scat¬ 
tering  angle  increases  the  excitation  light  incident 
on  the  detection  filter,  enhancing  spurious  dye 
measurements. 

Although  the  SCUFA  accurately  measures 
Rhodamine  WT  in  many  environments,  the  com¬ 
bined  effects  of  filters  and  scattering  angle  make 


15 

•© 

r 

- 1 - 

WET  Labs  ECO  Triple! 

i  10 

1 

-  Z)(dye) 

T  (turbidity. 

- 

q  5 

0 

. 

Fig.  2  Field  measurements  of  spurious  D  {black)  and  r 
{gray)  vs  time  with  no  dye  in  the  water  for  a  WET  Labs 
ECO  Triplet  and  b  Turner  SCUFA.  The  instruments  were 
deployed  midsurfzone,  0.5  m  above  the  seabed  in  1  m 
depth,  and  separated  1  m  in  the  alongshore  at  Scripps 
Beach,  CA,  with  small  spilling  waves.  Seaward  of  the  sur¬ 
fzone,  significant  wave  height  was  0.5  m  and  peak  period 
was  15  s.  Note  turbidity  units  (ntu,  ntu*)  are  similar  but  not 
equivalent  due  to  instrumental  differences  (Section  3.2) 
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the  SCUFA  dye  estimate  more  sensitive  to  sur- 
fzone  bubbles  and  sand  (spurious  dye  values 
sometimes  reached  80  ppb,  not  shown)  than  the 
ECO  Triplet.  Further  testing  to  determine  the  ef¬ 
fect  of  surfzone  bubbles  and  sand  on  ECO  Triplet 
dye  measurements  is  described  below  (Section  5). 


3  Surfzone  Dye  Fluorescence  Instrumentation 

In  situ  open-face-type  (frame-mounted,  fixed  lo¬ 
cation)  and  fiow-through-type  (mobile,  jet-ski- 
mounted,  Fig.  3)  fiuorometer  systems  were  used 
to  measure  Rhodamine  WT  dye  concentration 
(D)  in  the  surfzone.  Fixed-location  fiuorometers 
provide  much  higher  temporal  resolution  and 
longer  sampling  periods  than  previous  hand-filled 
bottle  methods,  but  extensive  spatial  coverage  is 
precluded  by  cost  and  logistics.  Increased  spatial 
coverage  is  provided  by  a  jet-ski-mounted  fiuo¬ 
rometer  system.  The  jet  ski  allows  fast  repeated 
cross-  and  alongshore  transects  though  the  sur¬ 
fzone  where  traditional  boats  and  submersibles 
cannot  go.  Simultaneous  turbidity  (r)  measure¬ 
ments  are  used  to  estimate  relative  bubble  and 
sand  interference  in  the  instrument  sample  vol¬ 
ume  and  to  correct  estimates  of  D.  These  in¬ 
struments  can  also  be  used  to  measure  surfzone 
chlorophyll-a,  but  the  methodology  differs  from 
that  of  dye  (Omand  et  al.  2009). 


Fig.  3  Schematic  of  jet  ski  mounted  WS  flow-through 
system  (not  to  scale) 


WET  Labs  ECO  Triplet  open-face-type  fiuo¬ 
rometers  (hereafter  ET)  were  used  to  measure 
dye  concentration  D  and  backscatter  turbidity  r 
(see  Section  2  for  specifications),  and  were  frame- 
mounted  at  fixed  locations  in  the  field.  The  ETs 
averaged  8  Hz  data  to  yield  sampling  rates  of 
0.89  Hz  in  the  laboratory  and  0.23  Hz  in  the  field. 

The  jet-ski-mounted,  fiow-through  fiuorome- 
ter/turbidity  sensor  system  (hereafter  WS)  con¬ 
sisted  of  a  WET  Labs  Wetstar  fiuorometer 
(470  nm  excitation  wavelength,  570  nm  fiuo- 
rescence  detection,  0.11-400-ppb  range)  and  a 
Turner  SCUFA  turbidity  sensor  (see  Section  2  for 
specifications)  with  fiow-through  cap.  Although 
the  SCUFA  is  capable  of  measuring  both  fiu- 
orescence  and  turbidity  (Section  2),  the  surf¬ 
zone  fluorescence  signal  was  noisy  (Fig.  2).  Here, 
the  SCUFA  (robust  enough  for  surfzone  con¬ 
ditions)  is  used  only  for  its  fiow-through  tur¬ 
bidity  measurement.  Water,  drawn  through  the 
jet-ski-mounted  intake  boom  from  20  cm  below 
the  surface  with  an  electric  pump,  passes  through 
a  debubbler,  the  SCUFA  turbidity  sensor,  and, 
finally,  the  Wetstar  fiuorometer  (Fig.  3).  The  de¬ 
bubbler  (a  small  200-ml  chamber)  allows  air  to 
escape  from  the  top  and  water  to  fiow  out  of  the 
bottom,  reducing  the  number  of  large  bubbles  en¬ 
tering  the  optical  instruments.  The  fiow-through 
system  hoses  and  debubbler  (Fig.  3)  smooth  sharp 
gradients  in  dye  occurring  over  times  less  than 
2.4  s,  and  delay  dye  measurements  relative  to 
GPS  positions  (Appendix).  Dye  fiuorescence  D, 
turbidity  r,  water  temperature  (measured  at  the 
end  of  the  intake  boom),  and  GPS  position  are 
sampled  at  5  Hz.  Data  are  both  logged  onboard 
and  transmitted  to  a  shore  station,  allowing  real¬ 
time  analysis  and  adaptive  sampling.  A  handlebar- 
mounted  screen  displays  real-time  data  and  a  local 
map  of  position,  allowing  repeated  sampling  of 
predetermined  transects. 

3.1  Fiuorometer  Calibrations 

Fiuorometers  were  calibrated  at  four  tempera¬ 
tures  (7-24  °C  range)  using  nine  known  concen¬ 
trations  of  Rhodamine  WT  (0-400-ppb  range)  in 
filtered  (to  remove  particles)  seawater.  Salinity, 
alkalinity,  and  pH  all  weakly  affect  dye  fiuores¬ 
cence  (Feuerstein  and  Selleck  1963;  Smart  and 
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Fig.  4  Laboratory  averaged  turbidity  r  vs  sand  concentra¬ 
tion  for  ET  {gray  circles)  and  WS  {black  triangles)  using 
Scripps  Beach  (La  Jolla,  CA)  sand.  For  the  open-face  ET, 
known  amounts  of  sand  were  completely  suspended  by 
vigorous  stirring  in  a  laboratory  tank,  and  turbidity  was 
averaged  over  2  min.  For  the  flow-through  WS,  sand-laden 
water  was  drawn  into  the  system,  and  the  average  dried 
sand  concentrations  were  measured  from  discharge  water 


Laidlaw  1977;  Stanbro  and  Pyrch  1979);  thus,  cal¬ 
ibrations  were  always  conducted  in  seawater.  Dye 
calibrations  were  found  to  be  within  15%  of  the 


factory  calibrations.  ET  calibrations  were  slightly 
nonlinear  above  100  ppb,  and  Wetstar  calibrations 
were  linear  up  to  400  ppb.  Calibrations  before 
and  after  the  field  deployment  were  similar.  Tem¬ 
perature  dependence  followed  Smart  and  Laidlaw 
(1977),  Dtq  =  Drexp[0.027(r  —  Tq)],  where  Dtq 
is  the  dye  fluorescence  normalized  to  a  reference 
temperature  To  (°C),  and  Dj  is  the  dye  fluores¬ 
cence  at  the  in  situ  temperature  T  (°C).  Factory 
turbidity  calibrations  were  used  for  both  the  ET 
and  WS. 

3.2  Turbidity  Measurements 

Turbidity  r  measured  with  the  ET  (ntu)  and  the 
SCUFA  (ntu*)  are  not  equivalent  because  the 
backscatter  wavelengths  and  angles  are  differ¬ 
ent  (Section  2).  In  laboratory  experiments  us¬ 
ing  Scripps  Beach  sand  suspended  in  seawater 
(Fig.  4),  SCUFA  and  ET  r  were  similar  for 
sand  concentrations  <  8  gL-i  (  <  30  ntu,  ntu*). 
However,  above  8  gL“^,  the  SCUFA  r  was  less 
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Fig.  5  Field  and  laboratory  results  without  dye  in  the  water 
for  {left  panels)  ET  and  {right  panels)  WS.  a,  c  Spurious 
D  vs  r  {solid  curves),  and  Diow  At  {dashed  curve)  from 
Eq.  1.  b,  d  Discrete  probability  density  functions  {PDF)  of 
r  corresponding  to  top  panels.  D  is  given  as  binned  means, 
with  vertical  bars  indicating  ±  one  STD.  ET  field  data 
{solid  black  curves)  combine  several  cross-shore  locations 


within  the  surfzone,  and  field  data  for  the  WS  combine 
multiple  cross-shore  transects.  Laboratory  data  are  given 
for  r  from  bubbles,  from  sand,  and  from  bubbles  and  sand 
combined  (see  legend).  Between  1,200  (lab)  and  28,000 
(field)  samples  are  used  to  construct  ET  curves,  and  1,180 
(field)  to  5,398  (lab)  samples  for  each  WS  curve 
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responsive  to  sand  than  the  ET.  The  SCUFA  r 
sensor  saturates  near  75  ntu*,  at  sand  concen¬ 
trations  >  20  gL“^  that  are  not  expected  above 
the  bottom  boundary  layer  in  the  surfzone.  Addi¬ 
tional  tests  (Fig.  5)  show  that  SCUFA  r  >  75  ntu* 
can  be  produced  by  bubbles  or  a  combination  of 
bubbles  and  sand,  but  not  by  sand  alone.  SCUFA 
T,  well  beyond  the  recommended  100  ntu*  upper 
limit  (Turner  Designs  SCUFA  manual,  revision 
2.1),  are  used  because  they  prove  useful  for  error 
estimation  and  dye  correction  in  the  WS  instru¬ 
ment  package  (Section  6).  The  r  response  to  bub¬ 
ble  void  fraction  was  not  measured. 


4  Instrument  Testing  Methods 

The  effects  of  surfzone-bubble  and  sand-induced 
turbidity  r  on  dye  measurements  D  were  ob¬ 
served  both  with  and  without  dye  in  the  water.  Ex¬ 
periments  without  dye  were  conducted  in  both  the 
field  (surfzone)  and  the  laboratory.  However,  con¬ 
trolled  dye  concentrations  could  not  be  produced 
in  the  field,  so  tests  with  known  (non-zero)  dye 
concentrations  were  performed  in  the  lab  where 
T  levels  were  varied  by  altering  the  amount  of 
suspended  bubbles  and  sand.  To  apply  laboratory- 
based  results  to  the  field,  it  is  assumed  that  natural 
surfzone  r  over  a  sand  bottom  is  primarily  caused 
by  bubbles  and  sand  (i.e.,  not  from  mud,  organic 
matter,  etc.). 

4.1  Field  Methods 

ET  and  WS  instruments  were  tested  over  sev¬ 
eral  days  in  the  surfzone  at  Huntington  Beach, 
California,  where  significant  wave  height  was 
about  1  m,  peak  periods  ranged  from  14  to  16  s, 
and  hourly  surfzone  averaged  alongshore  currents 
reached  0.5  ms“k  Huntington  Beach  sand  is  a 
tan,  medium-grained  quartz  typical  of  the  South¬ 
ern  California  coast.  Four  ETs  were  deployed  in 
various  cross-shore  locations  between  the  shore¬ 
line  and  4  m  mean  water  depth.  The  ETs  were 
mounted  nominally  0.5  m  above  the  bed  with  the 
sensor  facing  downward  at  30  degrees  from  verti¬ 
cal.  Collocated  temperature  measurements  were 
used  for  ET  dye  calibration  (Section  3.1).  The 
jet-ski-mounted  WS  was  driven  on  ~  200-m-long 


cross-shore  transects  from  seaward  of  the  surfzone 
towards  the  shore  (inbound  transect)  until  the  jet 
ski  turned  around  near  the  shoreline  (^  0.5  m 
water  depth)  and  returned  offshore  (outbound 
transect).  WS  field  data  were  corrected  for  delay 
time  (relative  to  GPS  positions)  within  the  flow¬ 
through  system  (Appendix). 

4.2  Laboratory  Methods 

Laboratory  experiments  with  and  without  dye  in 
the  water  were  conducted  for  both  the  open- 
faced  ET  and  the  flow-through  WS  system.  Black 
(to  reduce  reflected  light)  containers  were  filled 
with  15  1  of  filtered  seawater  and  five  dye  con¬ 
centrations  (0-238  ppb).  The  ET  was  held  in  the 
container  at  midwater  depth  (^  15  cm),  and  for 
the  WS,  the  intake  and  discharge  hoses  were 
mounted  near  the  bottom  of  the  container.  Bub¬ 
bles  (^0.1-5  mm  diameter),  roughly  represen¬ 
tative  of  those  injected  into  the  surfzone  by 
breaking  waves  (Deane  and  Stokes  1999),  were 
added  using  an  aquarium  aerator.  Scripps  Beach 
sand,  a  tan,  medium-grained  quartz,  was  added 
by  hand.  Each  system  was  perturbed  with  bub¬ 
bles;  then  with  sand;  and,  finally,  with  bubbles 
and  sand  combined.  Bubbles  and  sand  were  sus¬ 
pended  by  vigorous  stirring,  allowed  to  settle,  and 
then  the  process  was  repeated.  It  was  unknown 
if  the  bubble  and  sediment  size  distribution  of  a 
stirred  laboratory  container  closely  represented 
that  suspended  in  a  natural  surfzone.  However, 
the  amounts  of  bubbles  and  sand  in  the  laboratory 
were  varied  so  the  turbidity  range  spanned  that 
found  in  the  field. 


5  Results 

5.1  Field  and  Laboratory  Without  Dye 

Field  and  laboratory  data  {D  and  r)  without  dye 
in  the  water  were  collected  with  the  ET  (Fig.  5a, 
b)  and  the  WS  (Fig.  5c,  d)  instruments.  Spurious 
dye  concentrations  D  (when  no  dye  is  present) 
increase  with  turbidity  r  in  both  field  and  labo¬ 
ratory  measurements,  consistent  with  Smart  and 
Laidlaw  (1977).  High  r  events  are  infrequent  in 
the  surfzone  for  both  instruments  (Fig.  5b,  d). 
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resulting  in  low  mean  spurious  D  (0.25  ppb  for 
ET,  and  0.39  ppb  for  WS)  in  the  surfzone. 

The  standard  deviation  (STD)  of  spurious  dye 
measurements  (indicated  by  vertical  bars)  in¬ 
creases  with  T  for  the  ET  (Fig.  5a)  but  remains 
nearly  constant  for  the  WS  (Fig.  5c).  ET  labora¬ 
tory  D  results  are  similar  to  the  surfzone,  where 
combined  bubble/sand  results  are  the  closest  to 
the  field  (Fig.  5a).  WS  lab  results  are  also  similar 
to  the  surfzone  (Fig.  5c),  but  only  the  bubble- 
perturbed  test  produced  the  full  range  of  observed 
surfzone  turbidity  (because  of  the  strongly  nonlin¬ 
ear  T  response  to  sand.  Section  3.2). 

Spurious  dye  measurements  are  small;  how¬ 
ever,  the  integral  over  many  spurious  measure¬ 
ments  in  time  (ET)  or  space  (WS)  can  effect  dye 
concentration  statistics  (e.g.,  total  dye  mass,  dye 
patch/plume  width).  For  example,  if  the  WS  was 
driven  on  a  cross-shore  transect  without  dye  in 
the  water,  the  inclusion  of  spurious  dye  measure¬ 
ments  would  suggest  that  dye,  and  the  resulting 


ET 


T  (ntu) 

Fig.  6  Percent  change  in  measured  dye  %AD 
(Eq.  2)  vs  r  from  a,  d  bubbles;  b,  e  sand;  and 
c,  f  bubbles  and  sand  combined,  for  {left  panels)  ET  and 
{right  panels)  WS  with  four  dye  concentrations  in  the  lab 
(see  legend).  %AD  are  given  as  binned  means,  with  vertical 


plume/patch  width,  spanned  the  entire  transect. 
To  remove  spurious  dye  measurements  for  both 
instruments  in  the  surfzone,  a  lower  bound  for  dye 
detection  Diow  is  estimated  by  a  fit  to  D  >  99% 
of  the  spurious  surfzone  D  in  each  bin  vs  binned 
T.  The  disadvantage  of  this  approach  is  that  it 
sets  a  lower  limit  on  measurable  dye  concentra¬ 
tion  and  the  resultant  size  and  duration  that  a 
dye  patch/plume  can  be  observed.  However,  it  is 
suggested  that  D  below  Diow  be  set  equal  to  zero 
to  avoid  biasing  dye  statistics.  The  fit  (Fig.  5a,  c, 
dashed  black  line)  to  surfzone  field  data  without 
dye  for  each  instrument,  respectively,  is 

Aow  =  OfT  +  C  ,  (1) 

where  a  and  C  are  fit  constants  and  Diow  is  in 
ppb.  For  the  ET,  ot  =  0.02  ppb  •  ntu“^,  and  C  = 
0.3  ppb,  and  for  the  WS,  ot  =  0.01  ppb  •  ntu*“^,  and 
C  =  0.46  ppb. 


WS 


T  (ntu*) 


bars  indicating  ±  one  STD.  Between  550  and  860  ET  and 
3,600  and  5,400  WS  samples  are  used  to  construct  each 
curve.  Dashed  black  lines  in  c  and  d  show  fits  (Eq.  5) 
to  binned  mean  data 
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5.2  Laboratory  with  Dye 

Four  known  dye  concentrations  (Dq)  were  per¬ 
turbed  with  bubbles  (Fig.  6a,  d),  sand  (Fig.  6b,  e), 
and  a  combination  of  bubbles  and  sand  (Fig.  6c, 
f)  for  ET  (Fig.  6a-c)  and  WS  (Fig.  6d-f),  respec¬ 
tively.  Changes  in  measured  dye  concentration 
D  are  given  as  percent  change  %AD  from  the 
unperturbed  known  Dq 

%AD  =  100  (D  -  Do) /Do  .  (2) 

The  mean  %AD  is  almost  always  negative  (i.e., 
quenching)  and  is  similar  for  all  dye  concentra¬ 
tions  (curves  within  each  panel  of  Fig.  6  over¬ 
lay  each  other).  This  concentration-independent 
%  AD  is  consistent  with  Smart  and  Laidlaw  (1977) 
for  T -induced  quenching  of  Rhodamine  WT,  but 
contrasts  with  the  concentration-dependent  r- 
induced  quenching  of  chlorophyll-a  (Omand  et  al. 
2009).  The  STD  (vertical  bars)  about  the  mean 
%  AD  increases  with  r  for  the  ET,  and  to  a  lesser 
extent  for  the  WS.  The  relationship  between  ET 
%AD  and  r  is  similar  for  bubbles,  sand,  and 
bubbles  and  sand  combined  (i.e.,  curves  in  all  left 
panels  are  similar),  and  the  maximum  decrease 
(%AD  ^  -20%,  binned  mean  minus  one  STD) 
due  to  T  occurs  at  near  full-scale  turbidity  val¬ 
ues.  However,  bubbles  affect  the  WS  differently 
than  sand  (Fig.  6d,  e),  with  bubbles  producing  r 
up  to  300  ntu*  with  moderate  quenching  effects 
(binned  mean  %AD  minus  one  STD  as  large 
as  —25%),  and  sand  producing  r  up  to  75  ntu* 
with  large  quenching  effects  (binned  mean  %AD 
minus  one  STD  as  large  as  —40%).  The  large 
WS  sand-quenching  effect  at  moderate  r  values 
is  consistent  with  very  high  sand  concentrations 
(that  block  fluoresced  light)  required  to  approach 
the  75  ntu*  saturation  point  for  sand-induced  r 
(Fig.  4). 

6  Dye  Measurement  Corrections 

Dye  measurement  errors  are  reduced  by  dis¬ 
carding  dye  data  points  with  turbidity  r  above 
a  threshold,  and  correcting  the  remaining  dye 
for  T -induced  quenching.  ET  data  were  removed 
when  T  >  90  ntu  to  exclude  the  increased  scatter 


in  quenching  (Fig.  6c)  near  and  beyond  the  in¬ 
strument  T  saturation  point  (100  ntu).  WS  data 
were  removed  when  r  >  300  ntu*,  the  highest  lab- 
observed  T  (and,  thus,  the  limit  of  lab-based  r 
correction).  The  combination  of  discarding  high  r 
data  and  correcting  the  remaining  D  for  r -induced 
quenching  is  hereafter  called  correction. 

Surfzone  fluorometers  encounter  bubbles  and 
sand  suspended  by  breaking  waves  and  strong 
currents.  The  ET  turbidity  sensor  cannot  distin¬ 
guish  between  bubbles  and  sand.  However,  the 
effects  of  bubble-  and  sand-induced  r  on  ET  dye 
quenching  are  similar  (Fig.  6a,  b),  and  the  curve 
derived  from  the  combination  of  bubbles  and  sand 
(dashed  line  in  Fig.  6c)  is  used  for  ET  corrections. 
In  contrast,  bubbles-  and  sand-induced  r  produce 
very  different  dye-quenching  effects  in  the  WS 
(Fig.  6d,  e),  but,  as  with  the  ET,  the  instrument 
cannot  distinguish  between  r  from  bubbles  and 
T  from  sand.  In  the  present  application,  the  jet- 
ski-mounted  WS  system  samples  water  20  cm 
below  the  surface  where  average  sand  concentra¬ 
tions  are  expected  to  rarely  exceed  1  gL“^  (^-g-^ 
Beach  and  Sternberg  1992;  Yu  et  al.  1993;  Beach 
and  Sternberg  1996;  Ogston  and  Sternberg  2002). 
Even  extreme  sand  concentrations  (>  20  gL“^) 
only  result  in  r  up  to  75  ntu*  (Fig.  4)  with  the 
WS  T  sensor,  leaving  about  one  third  of  the  field 
data  with  higher  r  than  sand  alone  can  produce 
(Fig.  5d).  Therefore,  bubbles  are  assumed  to  be 
the  T  source,  and  WS  dye  measurement  correc¬ 
tions  are  made  using  the  bubble-quenching  curve 
(Fig.  6d). 

The  relationship  between  binned  mean  %AD 
(Eq.  2)  and  r  are  used  to  correct  for  r -induced 
dye  quenching.  The  fits  to  ET  bubbles  and  sand 
combined  (Fig.  6c,  dashed  black  line)  and  WS 
bubbles  only  (Fig.  6d,  dashed  black  line)  have  the 
form 

%AD  =  100^t.  (3) 

Substituting  Eq.  2  into  Eq.  3  yields 


where  the  fit  constant  is  p  =  1.2  x  10“^  ntu“^ 
(ET)  and  =  0.91  x  10“^  ntu*“^  (WS).  Equating 
Do  to  corrected  dye  concentration  Dc,  and  D  to 
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ET  WS 


Fig.  7  Laboratory  examples  of  a,  c  percent  change  in 
raw  %Ai)r  (dashed  black  curves)  and  corrected  %ADc 
(with  Eq.  5,  solid  black  curves)  dye  concentration  and  b, 
d  turbidity  r  vs  time,  for  (left  panels)  ET  and  (right  panels) 
WS.  Thick  dashed  gray  line  is  the  unperturbed  %AD  =  0. 


the  raw  dye  concentration  Eq.  4  yields  the 
equation  to  correct  for  r -induced  quenching 


In  laboratory  tests  with  known  Dq  =  10  ppb 
(Fig.  7),  the  magnitude  of  raw  %ADr  (Eq.  2) 
is  increased  by  r -induced  quenching  during  pe¬ 
riods  of  high  turbidity.  Corrected  dye  concen¬ 
trations  Dc  are  usually  more  accurate  than  raw 
Dc  (\%ADc\<\%ADr\).  Occasionally,  Z)r  is 
quenched  when  r  is  low  (e.g.,  at  ^  =  80  s,  Fig.  7a), 
and  significant  errors  (~  15%)  remain  in  cor¬ 
rected  Dc.  Discarding  ET  measurements  with  r  > 


For  the  ET  (WS)  seawater  with  a  known  Dq  =  70-ppb 
dye  concentration  was  perturbed  with  bubbles  and  sand 
(bubbles  only),  a  Light  gray  vertical  bands  indicate  times 
when  r  >  90  ntu.  At  t  ^  60  s,  the  ET  r  sensor  is  saturated 
at  100  ntu 


90  ntu  (light  gray  vertical  bars)  removed  the  large 
spike  at  t  =  60  s  (Fig.  7a)  when  the  r  sensor 
was  saturated.  Results  are  qualitatively  similar  for 
other  Do  (not  shown). 

Errors  in  raw  and  corrected  Dc  (discarding 
high  T  data  and  using  Eq.  5)  are  calculated  from 
laboratory  data  perturbed  with  a  combination  of 
bubbles  and  sand  (bubbles  only)  for  the  ET  (WS) 
with  four  dye  concentrations  (Table  1).  For  both 
raw  and  corrected  data,  percent  root-mean-square 
errors  %£rms  =  100  x  ((D^,c  -  / Dq  and 

mean  errors  %£:in  =  100  x  —  Dq)/Dq  are 
given  (Table  1),  where  {)  indicates  a  time  average. 
Raw  percent  error  magnitudes  are  independent  of 


Table  1  Percent  errors  (rms  %  firms,  and  mean  %fim)  in  the  laboratory  with  known  dye  concentration  Dq  (left  column) 
perturbed  with  bubbles  and  sand  (ET),  and  bubbles  (WS) 


Do  (ppb) 

ET 

WS 

%£^rms 

%fim 

%£^rms 

%fim 

Raw 

Corrected 

Raw 

Corrected 

Raw 

Corrected 

Raw 

Corrected 

35 

6.0 

5.0 

-3.9 

0.7 

12.1 

2.7 

-9.3 

-0.3 

70 

8.3 

4.9 

-6.7 

-1.6 

12.0 

1.8 

-8.9 

-0.1 

142 

8.5 

5.0 

-5.8 

-1.1 

12.4 

1.4 

-10.4 

-0.4 

238 

6.6 

3.3 

-4.7 

-1.0 

11.8 

2.6 

-8.9 

-1.8 

%  firms  and  %fini  are  given  for  raw  data,  and  data  corrected  by  removing  high  turbidity  (ET  r  >  90  ntu,  WS  r  >  300  ntu)  dye 
measurements  and  using  Eq.  5 
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flow-through  system  when  the  jet  ski  drives  over 
bubble-fllled  bores  (Fig.  10b).  On  inbound  tran¬ 
sects,  bubbles  are  reduced  by  driving  just  in  front 
of  the  shoreward  traveling  bore  (Fig.  10a).  WS- 
corrected  Z)c  on  sequential  inbound  and  outbound 
transects  have  closer  peak  values  than  D^. 

High-turbidity  events  associated  with  increased 
error  are  less  frequent  in  the  held  than  the  lab 
(Fig.  5b);  thus,  ET  dye  measurement  errors  in  the 
surfzone  are  expected  to  be  smaller  than  in  the 
laboratory  (Table  1).  In  contrast,  WS  held  errors 
may  be  increased  relative  to  laboratory  errors 
by  rare  high-concentration  sand  events  that  are 
unaccounted  for  in  lab  estimates. 


7  Summary 


Fig.  8  ET  a  raw  and  corrected  dye  concentration  (D^, 
Dc),  b  correction  magnitude  (Dc  —  D^),  and  c  r  vs  time 
in  the  surfzone.  The  instrument  was  mounted  65  m  from 
the  shoreline,  and  50  cm  above  the  bottom.  A  dye  patch 
(0.25  1  of  21%  by  weight  Rhodamine  WT)  was  released 
10  m  alongshore  from  the  ET  at  t  =  150s. a  Dc  (gray  curve) 
is  estimated  from  (black  curve)  and  c  corresponding  r 
time  series,  using  (Eq.  5).  Light  gray  vertical  bands  indicate 
times  when  r  >  90  ntu  (data  are  removed).  Light  gray 
hatching  indicates  times  when  <  D\q^  (Eq.  1),  and  Dr  c 
should  be  set  to  zero 

concentration  and  are,  generally,  <9%  (<  13%) 
for  ET  (WS).  For  both  ET  and  WS,  correction 
reduces  rms  (|%£rmsl  <  5%)  and  mean  (|%£ml  < 
2%)  errors.  Error  propagation  into  spatial  dye 
moments  is  nontrivial  and  depends  on  the  shape 
of  the  dye  distribution,  how  realizations  are  aver¬ 
aged,  and  the  assumed  noise  decorrelation  scale. 
However,  for  a  simple  Gaussian  example  with  5% 
rms  noise  and  a  noise  decorrelation  scale  equal  to 
one  STD  of  the  Gaussian  itself,  the  rms  error  in 
measured  variance  is  2%. 

Examples  of  corrections  to  surfzone  dye  mea¬ 
surements  are  shown  for  ET  (Fig.  8)  and  WS 
(Fig.  9).  ET  T -induced  quenching  corrections 
Dc  —  Dr  are  small  (<  4  ppb.  Fig.  8b)  compared 
to  dye  variability  (Fig.  8a).  WS  corrections  are 
similarly  small  (<  1  ppb  on  the  inbound  transect), 
but  they  reach  9  ppb  on  the  outbound  transect 
when  T  elevated  by  bubbly  water  is  drawn  into  the 


Open-face  (ET)  and  flow-through  (jet-ski- 
mounted  WS)  Rhodamine  WT  fluorometers  for 
in  situ  surfzone  sampling  were  tested.  Surfzone 


-80  -60  -40  -20  0 

cross-shore  coordinate  (m) 


Fig.  9  WS  a  D  and  b  r  vs  cross-shore  coordinate  (pos¬ 
itive  onshore,  with  shoreline  at  zero),  for  a  single  cross¬ 
shore  transect  pair,  575  m  alongshore  from  the  dye  source 
during  a  continuous  release  in  the  surfzone.  Inbound  and 
outbound  transect  raw  D  and  r,  and  corrected  D,  are 
shown  (see  legend).  Inbound  corrected  dye  curve  is  ver¬ 
tically  offset  by  -rl  ppb  for  visibility.  D  >  Diow,  and  root- 
mean-square  Diow  was  0.56  and  0.62  ppb  for  inbound  and 
outbound  transects,  respectively.  About  1  min  of  data  are 
shown 
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Fig.  10  a  Inbound  and 
b  outbound  cross-shore 
transects  through  the 
surfzone  with  the 
jet-ski-mounted  WS. 
a  Inbound  transect 
bubbles  are  minimized  by 
driving  just  in  front  of  a 
shoreward-traveling  bore, 
b  On  outbound  transects, 
bubbly  water  is  drawn 
into  the  flow-through  WS 
system  as  the  jet  ski 
drives  over  a  bore 


bubbles  and  sand  (measured  as  turbidity  r) 
interfere  with  fluorescent  dye  measurements  both 
raising  the  lower  bound  for  dye  detection  D\o^ 
(Eq.  1)  and  reducing  (quenching)  measured  dye 
fluorescence  D  (Fig.  6).  Laboratory  experiments 
with  known  dye  concentrations  Dq  in  seawater 
were  perturbed  with  bubbles,  sand,  and  bubbles 
and  sand  combined  to  explore  the  effect  of  r  on 
D.  The  reduction  in  measured  D  increases  with 
T,  and  the  observed  relationship  is  used  to  cor¬ 
rect  surfzone  dye  estimates.  Percent  root-mean- 
square  %£rms  and  mean  %Sm  errors  for  the 
ET  and  WS  are  reduced  to  <  5%  by  discarding 
high-T  data  and  using  Eq.  5.  In  situ  Rhodamine 
WT  measurements  are  feasible  in  the  surfzone, 
allowing  comprehensive  tracer  mixing  and  trans¬ 
port  experiments  in  this  heavily  used  but  poorly 
understood  region  of  the  ocean. 
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Appendix:  Flow-Through  Mixing 
and  Delay  Time 

Continuously  pumped  water  in  the  WS  flow¬ 
through  system  (Fig.  3)  is  mixed  (in  hoses  and 
the  debubbler),  smoothing  out  sharp  gradients  in 
dye  concentration.  In  addition,  water  takes  several 
seconds  (delay  time)  to  move  through  the  hoses 
and  debubbler  before  reaching  the  fluorometer. 
This  delay  time  must  be  accounted  for  to  match 
GPS  positions  with  dye  measurements.  Smoothing 
and  delay  time  are  specific  to  the  pump,  hoses 
(length,  diameter,  and  roughness),  and  debubbler 
(not  commercially  available)  used  in  the  WS  flow¬ 
through  system.  The  specific  smoothing  and  delay 
time  estimates  in  the  WS  system  are  presented 
here  to  provide  context  for  the  results  presented 


Fig.  11  Mean  (over  five  step  functions)  time  series  of  nor¬ 
malized  WS  dye  concentration  {black  curve)  in  response  to 
a  dye  step  function  input.  The  system  output  shows  a  10- 
90%  rise  time  (between  dashed  gray  lines)  of  2.4  s,  due  to 
mixing  in  the  flow-through  system  hoses  and  the  debubbler 
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inbound  —  —  -  -  inbound,  time  corrected 

outbound  -  -  -  -  outbound,  time  corrected 


Fig.  12  WS  dye  D  vs  cross-shore  coordinate  (positive 
onshore,  with  shoreline  at  zero),  for  single  inbound  (black) 
and  outbound  (gray)  cross-shore  transect  pairs,  a  50  m 
and  b  575  m  alongshore  from  the  dye  source  during  a 
continuous  release.  Time-corrected  inbound  and  outbound 
curves  (dashed)  are  similar,  suggesting  that  system  delay 
time  in  raw  data  (solid)  is  accounted  for 


above  and  to  establish  a  framework  for  character¬ 
izing  similar  systems  in  the  future. 

Smoothing  in  the  WS  flow-through  system  was 
estimated  with  laboratory  experiments  using  step 
functions  of  dye  (created  by  switching  the  intake 
between  water  with  zero  and  a  known  dye  con¬ 
centration  with  a  Y  valve,  Fig.  11).  The  2.4  s  10- 
90%  rise  time  gives  the  temporal  smoothing  scale. 
This  smoothing  in  time  gives  rise  to  smoothing 
in  space,  which  depends  on  jet  ski  speed  (i.e., 
2.4  s  multiplied  by  jet  ski  speed).  For  example, 
at  a  typical  jet  ski  speed  of  4  ms“^  (the  linear 
wave  speed  in  1.6  m  water  depth),  the  distance 
between  independent  samples  is  about  9.6  m.  The 
WS  flow-through  system,  essentially  a  low-pass 
filter,  results  in  smoothing  for  dye  concentration 
frequencies  greater  than  about  (2.4  s)“^  entering 
the  WS  system.  The  severity  of  the  smoothing 
increases  with  frequency  and  is  reduced  by  using 
slower  jet  ski  speeds  (thus  reducing  frequency). 
The  smoothing  magnitude  can  be  estimated  from 
the  assumed  true  dye  wavenumber  spectrum  and 
the  jet  ski  speed. 


The  WS  flow-through  system,  essentially  a 
low-pass  filter,  results  in  smoothing  for  dye  con¬ 
centration  frequencies  greater  than  about  (2.4  s)“^ 
entering  the  WS  system.  The  severity  of  the 
smoothing  increases  with  frequency  and  is  re¬ 
duced  by  using  slower  jet  ski  speeds  (thus  reduc¬ 
ing  frequency).  The  smoothing  magnitude  can  be 
estimated  from  the  assumed  true  dye  wavenumber 
spectrum  and  the  jet  ski  speed. 

Sequential  inbound-outbound  transect  pairs 
were  used  to  find  WS  system  delay  times  by 
minimizing  the  difference  between  the  inbound 
and  the  outbound  cross-shore  dye  concentrations 
(similar  to  Fig.  12).  From  64  transect  pairs,  a  delay 
of  8  s  was  the  most  common  value  (i.e.,  the  mode), 
with  a  mean  delay  of  8.1  s.  The  8-s  value  was  used 
to  correctly  time  the  WS  field  data.  In  examples 
of  WS  D  vs  cross-shore  distance,  sequential  in¬ 
bound  and  outbound  transects  are  closer  in  shape 
and  location  after  time  correction  (Fig.  12).  The 
addition  of  a  flow  rate  sensor  is  planned  to  more 
accurately  estimate  the  delay  time.  The  STD  of 
the  field  delay  times,  ±0.84  s,  results  in  spatial 
uncertainty  that  depends  on  jet  ski  speed.  For 
example,  the  spatial  uncertainty  is  approximately 
±3.3  m  at  4  m  s“^  speed. 
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